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ABSTRACT 


For  a subset  A of  Euclidean  n-space,  the  location  problem  of  finding  the 
point  x which  minimizes  the  maximum  distance  d(x,  a)  for  a in  A may  be  inter- 
preted as  finding  the  smallest  sphere  which  encloses  A.  Algorithms  have  pre- 
viously been  developed  for  the  situations  when  A is  finite  or  a polytope. 

Here  we  concentrate  on  interpretations  and  properties  of  the  problem,  particu- 
larly its  relation  to  other  problems:  cases  where  the  sphere  problem  is  dual 
to  that  of  finding  a shortest  vector  in  A,  where  there  is  a connection  with  a 
maximum  moment-of-inertia  problem,  etc.  Conjectures,  relating  the  minimum 
sphere  to  points  in  A which  define  the  diameter  of  A,  are  also  discussed. 


i 


1.  Introduction 


— 


For  a subset  S of  En,  the  location  problem  of  finding  the  point  x which 

minimizes  the  maximum  of  | |x  - a| | for  a e S,  may  be  interpreted  as  finding 

the  smallest  sphere  which  encloses  S.  Previously,  algorithms  have  been 

developed  for  solving  this  problem  when  S is  finite  or  a polytope  [4,  6].  In 

this  note  we  concentrate  on  interpretations  and  properties  of  the  problem; 

particularly  its  relation  to  other  problems.  For  example,  we  show  a case 

where  the  sphere  problem  is  dual  to  the  problem  of  finding  the  vector  of 

2 

minimum  norm  in  S . As  another  example  (in  E ) we  show  the  relation  to  a 
maximum  moment  of  inertia  problem.  Other  aspects  discussed  concern  conjectures 
relating  the  minimum  sphere  to  points  in  S which  define  the  diameter  of  S. 
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2 . Primal  and  Dual  Formulation 

The  minimum  sphere  for  a given  finite  set  of  points  a^  e En,  i = 1,  2,  . . . , m, 

may  be  determined  by  solution  of  the  mathematical  program 

(PI)  min  r 

(r,c) 

s.t.  r > | |c  - a.| | i = 1,  2,  ....  m 

with  variables  r and  c.  The  solution  of  (PI)  , (r*,  c*)  , gives  the  radius  and 
center  of  the  sphere  (the  existence  and  uniqueness  of  (r*,  c*) , for  any  bounded 
set  is  guaranteed  [1]).  Problem  (PI)  is  equivalent  to 
(P2)  min  s 

2 

s.t.  s > I |c  - ai| I i = 1,  2,  ....  m 

2 T 

where  s = r . By  introducing  the  variable  z = s - c c[9],  (P2)  becomes  the 

convex  quadratic  program 

( P 3)  min  z + c^c 

(z,c) 

T T 

s.t.  z + 2a.  c > a.  a.  i = 1,  2,  ....  m 
1=11 

The  dual  of  (P3)  is  the  concave  quadratic  program  [4]: 

(D3)  max  h(v)  = £ v^(a^a^)  ~ 

v 

s.t.  ^ v.  = 1 v.  > 0 i = 1 , 2 , . . . , m 

where  A has  columns  a.. 

i 

A short  note  by  Kuhn  [8]  shows  that  these  problems  (P3)  and  (D3)  are 

"hybrid"  programs,  i.e.,  they  are  linear  programs  with  a sum  of  squares  added 

to  the  objective  function.  (P3)  has  this  form  and  so  does  (D3)  if  we  define 

T 

w = Av  so  the  second  term  in  the  objective  is  w w. 

If  v*  solves  (D3)  and  (s*,  c*)  solves  (P2)  then 
s*  = h(v*) 

7 * 

and  c*  = ^ v.  a.. 
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r * 


Note  that  the  latter  equation  and  the  conditions  on  v in  (D3)  express  the 

intuitively  obvious  fact  that  c*  is  a convex  combination  of  the  a^.  The 

coefficients  expressing  the  optimal  combination  are  the  Kuhn-Tucker  multipliers 

which  solve  (D3) . By  a theorem  of  Caratheodory  [10],  the  convex  combination 

may  always  be  expressed  with  at  most  n + 1 positive  coefficients.  These  facts 

2 

are  easily  visualized  for  a.  e E . The  minimum  circle  will  either  be  defined 

y l 

by  the  two  most  distance  points,  or  it  will  pass  through  the  vertices  of  three 
points  which  define  a non  obtuse  triangle.  Of  course,  more  points  may  lie  on 
the  circle,  but  two  or  three  suffice  to  express  the  center  as  a convex  combination. 
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3 . Duality  Interpretations 

la  this  section  we  show  that  two  problems,  one  in  mechanical  physics  and 
the  other  a geometrical  "least  distance"  problem,  are  dual  to  the  minimum 
sphere  problem. 

The  first  of  these  might  be  called  a maximum  moment  of  inertia  problem. 

(I  am  grateful  to  my  colleague,  E.  J.  Muth , for  pointing  out  this  problem 
and  its  relation  to  the  sphere  problem.) 

It  is  a basic  principle  of  physics  that  a body  set  into  rotation  will  tend 
to  rotate  about  its  center  of  gravity.  This  happens,  for  example,  when  a 
frisbee  is  sent  spinning  through  the  air.  The  physical  principle  involved  is 
that  natural  forces  tend  to  minimize  the  rotational  energy  of  the  body  by 
choosing  the  center  of  gravity  as  the  axis  of  rotation. 

Consider  now  the  following  problem:  Suppose  we  wish  to  assign  masses  nr 
to  a finite  set  of  points  a^  on  a body  so  as  to  maximize  its  moment  of  inertia. 
We  assume  that  the  a^  are  connected  by  thin  rods  or  in  some  other  manner  in- 
volving negligible  mass.  For  example,  the  a^  may  be  points  on  a thin  lamina. 
Without  loss  of  generality,  we  may  assume  the  total  mass  to  be  distributed 
among  the  points  is  equal  to  one.  Then,  since  the  formula  for  moment  of 
inertia,  I,  is 

I = l m r 2 
L 1 i 

where  is  the  distance  from  the  point  a ^ to  the  axis  of  rotation,  we  may 
write  the  problem  as 

(PI)  max  l m^  f | a . — x I | 2 

mi 

s.t.  £ nu  = 1 m^i0  i = 1 , 2 , . . . , m 
where  x is  the  center  of  gravity.  But  owing  to  the  minimizing  feature  of  the 
center  of  gravity  we  may  write 
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(PI1)  max  min  £ m^jla^  " xl 


m.  x 
i 


s.t.  £ itk  = 1 , m^  ^ 0 i = 1 , 2 , . . . , m 

Thus  the  moment  of  inertia  problem,  (PI)  , is  a maximin  problem,  and  by  a well 
known  theorem  [10,  Corollary  37.3.2],  the  max  min  in  (PI1)  may  be  interchanged 
to  min  max.  This  yields 

2 

(Dll)  min  max  £ m^||a^  - x|| 


y m = 1 m.  > 0 i - 1,  2,  ...»  m. 
u i i = 


The  inner  maximization  is  easily  solved  by  setting  m.  = 1,  where  J j a . — x J | = 

2 

max^  | |a^  - x| | , thus  we  obtain 
(DI2)  min  max  | |a^  - x| 

Which  is  equivalent  to  (P2)  where  c = x and  c*  = x.  Thus  we  have  proven  that 

the  solution  of  the  problem  (PI)  is  obtained  by  constructing  the  minimum  circle 

* * 

for  the  given  a.  and  assigning  the  weights  m.  = v.  , where  v is  the  vector  of 
optimal  Kuhn-Tucker  multipliers. 

It  is  interesting  that  the  theorem  of  Caratheodory  quoted  earlier  states 
that,  for  this  problem,  the  maximum  number  of  points  with  positive  weight  need 
be  no  greater  than  three.  Also,  the  angular  kinetic  energy  and  the  radius  of 
gyration  of  a rotating  body  are  proportional  to  the  moment  of  inertia,  so 
maximizing  these  quantities  yields  the  same  result. 

3 

The  interpretation  given  above  does  not  extend  directly  for  a^  6 E because 
the  axis  of  rotation  is  a line  rather  than  a point. 

Our  second  problem  arises  in  the  minimization  of  nondif f erentiable  convex 
functions  [13].  An  important  subproblem  is  to  determine  the  vector,  d*,  of 


minimum  norm  in  a given  compact  convex  set  S,  i.e.,  solve  the  problem: 
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(PM)  min  ||d|i 

s.t.  d e S. 

If  we  make  the  assumption  that  S is  a polytope  with  extreme  points  a^ , 
i = 1 , 2 , ....  m,  then  we  have 

1 I I L viaill2 

V 


(PMl) 


Id  | I ^ = min 


s.t.  £ = 1 vi=>0  i = 1,  2,  ...»  m. 

If  we  further  assume  that  the  a.  have  the  same  Euclidean  norm , say  | |a.  j 

'L  i 

then  from  (PMl)  we  obtain 


1 - ||d*||2  - 1 - »in  III  v,a,||2 


= 1 + max  - ! I T v . a . | 

11 L i l 1 


= max 


1 ~ I II  ^aj 


= max  l vi(aiTai)  - ||]>  v^JI2 
s.t.  I v£  = 1,  v.  > 0. 

Comparison  with  (D3)  gives  the  (Pythagorean)  dual  relation 
||d*||2  + s*  = 1 

where  s*  is  the  radius  squared  of  the  minimum  sphere  covering  the  a.  and  d*  = c* . 
See  Figure  1 for  a sketch.  See  [ 7 ] for  computational  experience  where  such 
normalizing  of  the  a^  has  proven  helpful  in  nondif ferentiable  optimization. 

It  would  be  interesting  to  know  whether  there  is  a geometrical  interpretation 
for  the  case  when  not  all  | |a^| | = 1. 
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As  an  aside  on  problem  (PM1),  we  note  that  it  may  be  written  as 


T T 

(PM2)  max  [ - d d + 2 min  a.  d]. 
d i 1 

This  formulation  calls  for  the  unconstrained  maximization  of  a strictly 
concave  function,  but  as  far  as  we  know,  it  has  not  been  exploited  in  any 
solution  procedure.  To  obtain  (PM2)  from  (PMJL),  one  has  only  to  regard 


4.  Solution  Methods 


Since  both  (P3)  and  (D3)  are  quadratic  programs,  either  could  be  solved 
by  one  of  the  well  known  algorithms  designed  for  such  problems.  Care  should 
be  taken,  however,  because  both  objectives  are  semide finite.  Thus  methods 
such  as  the  Theil-Van  de  Panne  procedure  [12 j which  require  a definite 
quadratic  objective  are  not  directly  applicable,  but  others,  such  as  Dantzig's 
procedure  [ 2 ] are.  If  m is  small,  then  (P3)  or  (D3)  could  both  be  solved 
with  about  the  same  effort.  If,  however,  m is  very  large,  (D3)  is  more  attrac- 
tive because  there  is  only  one  constraint.  In  [ 4 ] Elzinga  and  Hearn  give 
a decomposable  solution  procedure  for  (D3)  which  uses  the  Dantzig  algorithm 
and  for  which  the  size  of  the  tableau  is  independent  of  m.  The  basis  for  the 
decomposition  is  Caratheodory ' s theorem.  A similar  algorithm  for  the  case 
when  S is  a polytope  is  described  in  [ 6 ]. 

For  the  problem  in  the  plane  one  could  also  employ  quadratic  programming, 
but  the  following  geometrical  procedure  quoted  from  [ 5 ] is  more  efficient: 

"1.  Pick  any  two  given  points  and  go  to  step  2. 

2.  Let  the  two  points  define  the  diameter  of  a circle.  If  this  circle 
covers  all  points,  stop.  Otherwise  choose  some  point  outside  the 
circle  and  the  two  defining  points  and  go  to  step  3. 

3.  If  the  three  points  define  a right  triang]  or  an  obtuse  triangle, 
drop  the  point  at  the  angle  ^90°  and  go  to  ■* cp  2 with  the  remaining 
two  points.  Otherwise  if  the  triangle  has  strictly  acute  angles,  go 
to  step  4. 

4.  If  the  circle  defined  by  the  three  points  covers  all  points,  stop. 
Otherwise  choose  some  outside  point,  call  it  a,  and  label  as  a a 
point  from  among  the  three  defining  points  that  is  farthest  from  a 
Extend  the  diameter  of  the  current  circle  through  point  a to  divide 
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the  plane  into  two  half-planes.  Label  as  3 the  defining  point  in  the 
same  half-plane  with  o and  as  y the  remaining  point.  With  points  a,  y, 
and  a,  go  to  step  3. 


,! 


! 

I 


El 


In  step  3 if  the  three  points  are  on  a straight  line,  we  con- 
sider this  a 'triangle'  with  a 180°  angle  at  the  intermediate  point. 
Also  note  that  by  considering  a right  triangle  as  a two-point  case 
we  assume  that  in  step  4 the  point  3 is  strictly  contained  in  the 
half-plane  with  a.  In  practice  it  seems  reasonable  to  choose  the 
farthest  point  outside  the  current  circle  rather  than  just  some 
point  as  indicated  in  steps  2 and  3. 

Convergence  of  this  algorithm  is  based  on  showing  that  the 
circles  created  in  each  iteration  are  monotone  increasing  in  radius. 
Then  since  there  are  only  finitely  many  two-point  and  three-point 
circles,  the  process  is  finite." 

We  quote  the  algorithm  here  (but  not  the  convergence  proof)  to  emphasize 
(a)  that  it  is  initiated  with  an  arbitrary  pair  of  points,  and  (b)  that  once  a 
point  is  covered,  it  need  not  remain  covered  (see  Figure  2).  All  the  pro- 
cedure guarantees  is  that  successive  circles  strictly  increase  in  radius. 


Point  a^  uncovered  by  third 
circle. 

FIGURE  2 
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The  efficiency  of  this  geometrical  algorithm,  as  measured  by  empirical 

testing,  appears  to  be  good  for  m < 100.  (See  Figures  3 and  4.)  For  large  m, 

2 

however,  the  time  is  probably  of  order  m ; in  fact,  this  is  suggested  by 
Figure  4. 

Whether  a more  efficient  method  can  be  devised  is  an  open  question. 

Certainly,  an  algorithm  which  kept  all  prior  points  covered  would  be  of  order 

m*,  but  would  require  that  each  circle  be  efficiently  constructed  from  the 

previous  one.  In  other  words,  the  critical  question  would  be  how  to  construct 

the  minimum  circle  for  k + 1 points,  knowing  the  minimum  circle  for  k points. 

The  answer  is  trivial  for  k + 1 = 2 or  3,  but  for  k + 1 = 4,  there  are 
A A 

(3)  + ^ = 10  possible  circles  to  consider,  and  for  larger  k the  subproblem 
seems  no  easier  than  the  original. 


* Shamos  and  Hoey  [11]  have  conjectured  that  an  algorithm  of  order  m log^  m 
is  possible. 
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5 . The  Minimum  Circle  and  Far  Points 

2 

For  a finite  point  set  S = {a^,  a^9  a^},  ai  e E , we  define  as  far 

[Joints  any  pair  (a^,  a^)  for  which  | | a^  - a^|  | = max  | | a^  - a.j|  = diameter  of  S. 

i»j  J 

It  is  natural  to  conjecture  that  the  minimum  circle  passes  through  at 
least  one  far  point.  Figure  5 represents  a counterexample.  The  far  points  are 
and  a^  because  the  two  arcs.  Arc  1 and  Arc  2 enclose  all  a^.  Clearly,  however 
the  minimum  circle  is  defined  by  a^,  a^,  and  a^.  (The  distances  from  a^  to 
a^  and  a^  are  less  than  the  diameter  of  the  circle.) 


Minimum  Circle  defined  by 
{a3’  V a5}* 

FIGURE  5 

The  same  figure  enables  us  to  determine  the  probability  that  a circle 
defined  by  two  points  will  enclose  a random  third  point. 
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Theorem  Let  and  a ^ be  fixed  points  and  assume  a random  third  point  is  no 
farther  from  a^  and  a ^ than  | |a^  - a^  | | . Then  the  probability  P,  that  the 
circle  defined  by  and  a ^ covers  the  third  point  is 


P = 


3ti 


.64. 


8 - 6/  3 

Proof:  Consider  Figure  6 where  a^  and  a^  (distance  d apart)  represent  the  two 
fixed  points  so  that  the  third  is  enclosed  in  the  area  bordered  by  Arc  1 and 
Arc  2.  The  probability  of  it  not  being  exterior  to  the  circle  defined  by  a^ 
and  a^  is  simply  the  ratio  of  the  area  of  the  circle  to  the  area  between  the 
arcs. 

Point  a^  is  placed  at  the  origin  and  the  arcs  each  have  a radius  of  d. 

We  determine  the  area  between  the  arcs  by  integration  over  one-quarter  of  the 
area  (shaded  portion): 


9,2  t •r3-  4.  • -1  «r3~ 

= 2d  ( - + sin  2 - 2 > 


= 2d  ( 


2 , it  /IT. 


4 * ‘ 


Since  the  area  of  the  circle  is  tt  , the  ratio  is 


TT 


3/2  TT 


3tt 


q , tt  / 3 . 4tt  - 3/  3 8tt  - 6/  3 

8 ( 3 - ^ ) 


.64 


f(x) 


/7T7  - 4 


Arc  1 


Arc  2 


(d,  0) 


FIGURE  6 


This  result  suggests  further  problems  of  the  same  sort.  For  example,  what 

is  the  probability  that  three  points  uniformly  distributed  over  some  convex  set 
2 

in  R will  form  a triangle  with  an  obtuse  angle  (so  that  the  minimum  circle  is 
defined  by  the  two  far  points)?  When  the  set  is  the  unit  square,  Monte  Carlo 
simulation  shows  that  the  answer  is  approximately  .70.  To  our  knowledge,  the 
question  posed  has  not  been  answered  analytically.  For  the  interested  researcher, 
the  monograph  by  Kendall  and  Moran  [14]  is  a starting  point. 
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